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ABSTRACT 

We report the discovery of a new binary pulsar, PSR J1829+2456, found during a 
mid-latitude drift-scan survey with the Arecibo telescope. Our initial timing obser- 
vations show the 41-ms pulsar to be in a 28-hr, slightly eccentric, binary orbit. The 
advance of periastron uj = 0.28 ± 0.01 deg yr~^ is derived from our timing observa- 
tions spanning 200 days. Assuming that the advance of periastron is purely relativistic 
and a reasonable range of neutron star masses for PSR J1829+2456 we constrain the 
companion mass to be between 1.22 M© and 1.38 M©, making it likely to be another 
neutron star. We also place a firm upper limit on the pulsar mass of 1.38 Mq. The 
expected coalescence time due to gravitational-wave emission is long (~ 60 Gyr) and 
this system will not significantly impact upon calculations of merger rates that are 
relevant to upcoming instruments such as LIGO. 
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1 INTRODUCTION 

The first binary pulsar B1913-I-16 was discovered by Hulse & 
Taylor (1975). This double neutron star (DNS) system with 
its 7.75-hr orbital period and large eccentricity (e — 0.6) 
has since become a wonderful laboratory for testing gen- 
eral relativity in the strong-field regime (Taylor & Weis- 
berg 1989). Perhaps most importantly, it has provided the 
first evidence for the ex istence of gravitational radiation 
jTavlor fc WeisberdlTgsi) . DNS binaries start life as binary 
systems with main-sequence stars of mass > 6 Mq . Eventu- 
ally the more massive of the two stars undergoes a supernova 
explosion, leaving a neutron star, sometimes seen as a pul- 
sar. As the less-massive star evolves it increases in size until 
it overfills its Roche lobe. At this point, matter starts to ac- 
crete onto the pulsar causing it to spin up to periods as short 
as a few milliseconds (Alpar et al. 1982), a process known 
as recycling. In most cases the outer layers of the compan- 
ion are blown away after accretion, exposing the core of the 
companion and producing a white dwarf-millisecond pulsar 
binary. For some recycled systems, however, the companion 
is massive enough to explode as a supernova. In most cases, 
this violent event will disrupt the binary system. The DNS 
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binaries are those systems fortunate enough to survive. F or 
further details see iBhattacharva fc van den HeuveiJl99J) . 

Our understanding of the DNS binary population is 
currently hampered by small-number statistics. Large-scale 
pulsar surveys are being carried out by a number of groups 
in order to improve this situation by increasing the sam- 
ple of objects. In this Letter we report on the results of a 
drift-scan survey using the Arecibo telescope. The survey 
was conducted at 430 MHz and covered mostly intermedi- 
ate Galactic latitudes < 60°) to optimise the likelihood 
of finding millisecond and binary pulsars (see e.g. Cordes fc 
Chernoff 1997). The known pulsars detected are described 
along with preliminary parameters for a new 41-ms pulsar 
J1829-I-2456 which is likely to be a DNS binary. We com- 
pare the properties of this new system with the known DNS 
and other relativistic binaries. Finally, we outline the future 
observational prospects for this system. 



2 OBSERVATIONS AND ANALYSIS 

In July and August of 1999 the Arecibo radio telescope 
was used to take the data presented here. The receiver was 
parked and observations were made as the sky drifted over- 
head. The newly commissioned 430-MHz receiver system 
in the Gregorian dome was used in combination with the 
Penn State Pulsar Machine (PSPM), a 128-channel ana- 
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logue filterbank spectrometer which samples the incoming 
voltages fro m the telescope every 80 over a bandwidth 
of 7.68 MHz ^^^^M)- The Gregorian dome 430-MHz 
receiver has a gain of 11 K Jy~^ and a system tempera- 
ture of 45 K. At these latitudes a source passes through the 
beam in ~ 40 s. This is well matched to a 2^^-sample time 
series collected with the PSPM (42 s). Our survey sparsely 
sampled a total of ~ 230 deg^ within the right ascension 
range: 15h < a < 20h for two declination ranges defined by 
8° < 5 < 12° and 2A° < S < 28°. The instrumental setup 
was identical to that described by Lorimer et al. (2004), in 
which the system sensitivity was estimated to be greater 
than or about 0.5 mjy for long period pulsars away from 
the Galactic plane and 3 mJy for millisecond pulsars. Full 
details of the search sensitivity will be given elsewhere. 

The observations were analysed at the Jodrell Bank Ob- 
servatory using a 182-processor Beowulf cluster (COBRA 
and an identical procedure to that descr ibed by Lorim er et 
al. (2004). A suite of data analysis tools iLorimeill200 jl was 
used in combination with scripts to keep each node con- 
tinuously processing from the pool of data. Each drift ob- 
servation was split into 42-s beams, with each beam over- 
lapping 50% of the previous beam. These beams were then 
dedispersed using trial dispersion measures (DMs) between 
and 491.2 pc cm~^. The resulting time series were then 
Fourier transformed and searched for periodic signals. Inco- 
herent summing of the first 2, 4, 8 and 16 harmonics was 
used to increase sensitivity to narrow pulse profiles. Any re- 
sulting candidates with a signal to noise ratio (S/N) greater 
than 8 were folded for visual inspection. It took ~ 13,000 
CPU hours to complete the analysis of the ~ 14,200 over- 
lapping beams. 



3 RESULTS AND FOLLOW-UP 
OBSERVATIONS 

In the ~ 230 deg^ of sky covered by the survey, 4 known pul- 
sars were detected; PSR J1649-I-2533 (period, P = 1.02 s; 
DM = 34 pc cm"^) with S/N ~ 80, PSR J1652+2651 
(P = 0.92 s; DM = 40.8 pc cm~^) with S/N ~ 73, 
PSR J1532-h2745 (P = 1.124 s; DM = 14.7 pc cm"^) 
with S/N ~ 38 and PSR J1543-H0929 (P = 0.75 s; 
DM = 35.2 pc cm"^') with S/N ~ 85. The rate of pulsar 
detections in the data is consistent with the rate in previ- 
ous drift-scan data taken using the line-feed in combination 
with the PSPM (McLaughlin et al. 2003). 

One promising pulsar candidate was found with P = 
41.02 ms, DM ~ 13.9 pc cm"^ and S/N ~ 10. The candidate 
was reobserved and confirmed as a new pulsar, J18294-2456, 
on June 3 2003 using the same observing system with a 10- 
min integration. Further observations showed that the pul- 
sar scintillates strongly, with the flux density ranging from 
0.07 mJy to 0.7 mJy. In Fig. Q we present an integrated 
pulse profile. The confirmed approximate po sition put the 
pulsar in the line of sight of the Gould belt llGehrels et alJ 
bOOdMGrenieJbOOdl . a dense region of atomic and molec- 
ular gas containing young stars and an ideal environment 
for producing strong scintillation. The DM from the search 
observation and the pulsar's Galactic coor dinates (I ^ 53°; 
b ~ 15.6°) imply a distance of ~ 1.2 kpc iCordes fc Laziol 
|20o3). The true distance may in fact be much lower due to 
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Figure 1. Integrated pulse profile, formed from 4 hours of data, 
for PSR J1829-(-2456 at 430-MHz. 



the unmodelled excess electron density in the Gould belt. 
The confirmation period was 40.99 ms, significantly shorter 
than the discovery period. This indicated that the pulsar 
was most likely a member of a binary system. 

In order to determine the orbital parameters of 
PSR J1829-I-2456, observations were made using the same 
observing system over the following 6 months. These data 
were dedispersed and searched for the period of the pul- 
sar. Long observations when the pulsar appeared strongly 
were split and searched independently to provide more data 
points. The topocentric periods and observation times were 
then converted to barycentric periods and times using the 
TEMPO^ software package. These barycentric periods and 
times were then fitted to possible binary orbits, resulting in 
the curve shown in Fig. |5| which represents the pulse pe- 
riod evolution for an assumed orbital period P(, ~ 28 hr, 
projected semi-major axis x = ai sin i ~ 7.2 It-s and eccen- 
tricity e ~ 0.14. 

Using this preliminary ephemeris, the data were folded 
modulo the predicted period and cross correlated with a high 
S/N template profile to obtain an accurate time of arrival 
(TOA) for each observation. The TOAs were then analysed 
with TEMPO to produce a preliminary phase-connected so- 
lution (where every rotation of the pulsar is accounted for) 
over the time span of the observations. This resulted in the 
ephemeris given in Table Q 



4 DISCUSSION 

4.1 Mass determination 

The orbital parameters in Table can be used to constrain 
the mass of the binary system. The Keplerian mass function 
relates the orbital period. Pi,, and the projected semi-major 
axis, X, to the masses of the binary components. For this 
system we calculate a mass function 



/(mi,m2) = 



An'^x^ _ (m2 sini)^ 



(mi + 1712)'' 



= 0.294 M, 



01 



(1) 
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Figure 2. The orbital fit (line) to the barycentric periods and 
times (points) folded at the binary period of the system. The error 
bars are smaller than the points at this scale. 



where Tq = GMqc~^ = 4.925490947/iS, i is the inchnation 
between the plane of the orbit and the hne of sight, x is in 
Hght-seconds, Pb is in seconds and the pulsar and companion 
masses mi and m2 are in Solar masses. Neutron stars are 
ob served to have a narrow range o f masses. Given the data 
in iThorsett fc Chakrabartvl il999l) , and the low measured 
mass of PSR J0737-3039B (Lyne et al. 2004), we expect the 
pulsar mass to lie between 1.25 M0 and 1.47 Mq. For an 
orbit viewed edge-on (i = 90°), and a minimum neutron star 
mass for PSR J1829-I-2456 of mi = 1.25 M©, we calculate 
a minimum companion mass m2 — 1.22 M©. As shown in 
the following equation, the measurement of the advance of 
periastron, iu, if assumed to be purely relativistic, allows a 
measurement of the sum of the masses: 



(2) 



When the measured uj is used in combination with the con- 
straint that sini < 1, we obtain an upper limit on the 
maximum pulsar mass mi < 1.38 M©. Similarly, assuming 
mi > 1.25 Mq, the maximum value of Cj implies i > 66°. 
These constraints are shown in a mass-mass diagram in 
Fig. 121 Although we cannot currently rule out a massive 
white dwarf or main-sequence star, given the similar mass 
functions, spin and orbital parameters of J 1829+2456 to 
other DNS binaries (see Table it seems most likely that 
the companion is another neutron star. 



4.2 A search for the companion 

The accuracy of the position of PSR J1829-I-2456 inferred 
through radio timing is sufficient to allow for a search for 
an optical companion. No optical companion is present in 
the uncalibrated plates of the Digitised Sky Survey. Prepa- 
rations are underway for a deep optical search for the com- 
panion star. 

Given the possibility that the companion of J1829-I-2456 
is a neutron star, a search was made of the available radio 
data for a periodic signal which could be coming from the 
companion if it were active as a radio pulsar. This search was 
given greater significance following the recent detection of 
2.8-s pulsations from the companion of PSR J0737— 3039A 



1.4 



1.2 



\^cj niax 




1 1 1 1 






^^-\= 66°- 


\cj min\ 








^^-^1 = 90°" 


1 




1 Y. , N, 



CO 

a 

o 

'3 

e 

o 
u 



0.8 

1 1.5 

Pulsar mass (Mg) 

Figure 3. The constraints on the masses of the pulsar and com- 
panion. The vertical lines are the mass limits for neutron stars 
described in Thorsett & Chakrabarty (1999), and the low mea- 
sured mass of PSR J0737-3039B (Lyne et al. 2004). The lower 
companion mass limit is given by the mass function assuming a 
maximum inclination of 90°. The advance of periastron, tl), if as- 
sumed to be purely relativistic, provides the constraint for the 
maximum companion mass and maximum pulsar mass. 



(Lyne et al. 2004). Initial searches of the data (at the DM of 
PSR J 18294-2456) looking for other periodic signals showed 
that, if the companion was an active radio pulsar, it was too 
weak to be easily detectable. 

To carry out a more sensitive search for the compan- 
ion, we made use of the fact that the orbital parameters 
of the system are known and that the pulsations from the 
companion star would be Doppler shifted in the opposite 
sense to PSR J 1829-1-2456 due to its orbital motion about 
their common centre of mass. To correct for this effect, we 
applied the first-order Doppler formula to calculate the ef- 
fective sampling interval in the rest-frame of the companion: 



^companion — ^samp 



(3) 



where tsamp is the sampling interval at the observatory and 
v{€) is the effective line-of-sight radial velocity of the com- 
panion as a function of time, t. We used TEMPO to calcu- 
late «(t) to account for the effects of orbital motion in the 
binary system (using the orbital parameters given in Ta- 
bleQ assuming a pulsar mass of 1.35 M0 and changing only 
the longitude of periastron by 180°) as well as contributions 
from the Earth's rotation and motion about the Sun. All 
available data were dedispersed and resampled in this way 
before being passed through the standard search analysis. 
No significant candidates were found down to a S/N of 5. 
This corresponds to a 430-MHz flux limit of ~ 0.1 mjy. 

While the lack of detection of the companion suggests 
that it is too weak to be seen as a radio pulsar, or un- 
favourably beamed, it should be noted that the companion 
to PSR J0737— 3039A is only clearly detectable over certain 
phases of the orbit (Lyne et al. 2004). If any companion to 
PSR J 1829-1-2456 was behaving in a similar fashion it may 
not be visible in the current data set. Further, more sensi- 
tive, searches will be carried out as more data are taken. 



4 D. J. Champion et al. 



Table 1. Measured, derived, and expected parameters for 
PSR J1829+2456 from 29 TOAs from MJDs 52786 to 52986. 



Measured Parameters 



Right ascension (J2000) (h:m;s) 


18:29:34.6 (1) 


Decimation (J2000) ( : : j 


24:56:19 (2) 


Period (ms) 


41.00982358 (1) 


Epoch of period (MJD) 


52887 


Projected semi-major axis (s) 


7.2360 (1) 


Dispersion measure (pc cm~'^) 


13.9 (5) 


Binary period (days) 


1.176028 (1) 


Eccentricity 


0.13914 (4) 


Longitude of periastron (deg) 


229.94 (1) 


Epoch of periastron (MJD) 


52848.57977 (3) 


430-MHz flux density (mjy) 


0.3 (1) 


Pulse width at 50% of peak wrM (ms) 


1.07 (16) 


Pulse width at 10% of peak uiio (ms) 


2.07 (16) 


Rate of advance of periastron (deg yr~^) 


0.28 (1) 


Mass function (Mq) 


0.29413 (1) 


RMS residual to fit (^s) 


19 


Derived Parameters 


Galactic longitude (J2000) (deg) 


53.343 (1) 


Galactic latitude (J2000) (deg) 


15.612 (1) 


Distance (kpc)" 


~ 1.2 


Mean orbital speed (km s~^) 


134 


Total system mass (Mq) 


2.5 (2) 


Minimum companion mass (Mq) 


1.22 


Maximum pulsar mass (Mq) 


1.38 


Expected Relativistic Parameters 


Orbital period derivative (xlO"^'^) 


-0.02 


Coalescence time (Gyr) 


60 


Geodetic precession rate (deg yr~^) 


0.075 


Relativistic time-dilation 




and gravitational redshift 7 (ms) 


1.3 


a: Distance inferred from DM ICordes Lazi' 


J2002h. 



The numbers in parentheses are the l-tr uncertainties in the least 
significant digit quoted. 

4.3 Expected relativistic parameters 

Using the constrained range of masses for the pulsar and 
companion the expected orbital period derivative and period 
for geodetic precession due to misalignment of the spin axis 
of the pulsar with the orb ital angular momentum vector 
jBarker fc O'ConnellllToT^) can be determined. Assuming 
a combined mass of 2.53 Mq, the expected values for Pi, 
and the expected timescale for geodetic precession are given 
in Table The expected 7 should be measurable after 18 
months of timing data. The expected Pb is too small to be 
measurable without a long-term (5-10 yr) timing campaign. 
The expected timescale for geodetic precession is ~ 4200 yrs, 
too long to be measurable through pulse profile or timing 
variations. 

4.4 Comparison with other relativistic binary 
systems 

If our future measurements show PSR J1829-I-2456 to be a 
double neutron star system it will be the seventh such sys- 
tem to be discovered. Table |5| contains a list of the most rel- 
ativistic binary pulsar systems and their orbital parameters. 



We also list the coalescence time tqw due to the emission of 
gravitational radiation of each system. Following Lorimer 
(2001)'^, we can approximate the detailed calculations of 
TGw by Peters (1964) via 

rGw-10yr(-j — j (^^^ j ) / , 

where the reduced mass /i = mim2/(mi + 1112) ■ For any 
reasonable range of mi and m2 discussed above, we find 
TGW ~ 60 Gyr. 

The systems listed in Table |5| span a wide range of 
orbital parameters but can be broadly split into systems 
that will and will not coalesce within a Hubble time, 
with PSR J 1829-1-2456 lying toward the edge of the non- 
coalescing group. The relatively long coalescence time for 
J1829-I-2456 means that it will not affect the DNS merger 
rate calculations (e.g. Kalogera et al. 2004) for gravitational 
wave detectors such as LIGO. 



5 FUTURE STUDY 

In an Arecibo drift-scan survey of 230 deg'^ we have dis- 
covered a new relativistic binary system, PSR J1829-I-2456. 
Our measurements so far suggest that the companion is most 
likely another neutron star. The long predicted gravitational- 
wave coalescence time means that this system has no effect 
on the calculations of merger rates required for instruments 
such as LIGO. A dedicated Arecibo timing campaign will 
allow more accurate measurements of the position of this 
system, the spin-down rate of the pulsar and the relativistic 
advance of periastron. Continued timing will provide ad- 
ditional valuable information: the system's proper motion, 
the relativistic time-dilation and gravitational redshift pa- 
rameter 7, and possibly a measurement of Shapiro delay. 
An optical search for the companion is planned to confirm 
whether it is a massive white dwarf or, as seems most likely, 
another neutron star. Studies of scintillation effects are also 
planned which should allow the inclination angle of the orbit 
to be determined. 
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2 We reproduce the formula here to correct a typographic error 
in the original citation. 
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of various 


eccentric binary systems. 
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Pb 


ai sin i 
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/(m) 


mi ^ 


- m2 


rl^ References 




(ms) 


(days) 


(It-s) 




(degyr-i) (xlO^^^) 


(Mq) 


(M 


0) 


(Gyr) 




Double neutron star binaries 


B1913+16 


59.03 


0.323 


2.34 


0.617 


4.227 -2.428 


0.13 


2.83 




0.31 


1 


B1534+12 


37.90 


0.421 


3.73 


0.274 


1.756 -0.138 


0.31 


2.75 




2.69 


2 


B2127+11C 


30.53 


0.335 


2.52 


0.681 


4.457 -3.937 


0.15 


2.71 




0.22 


3 


J1518+4904 


40.93 


8.634 


20.04 


0.249 


0.011 


0.12 


2.62 




9600 


4 


J1811-1736 


104.18 


18.779 


34.78 


0.828 


0.009 <30 


0.13 


2.6 




1700 


5 


J0737-3039A 


22.70 


0.102 


1.42 


0.088 


16.88 -1.24* 


0.29 


2.58 




0.087 


6 


White dwarf binaries 


B2303+46 


1066.37 


12.34 


32.69 


0.66 


0.010 


0.25 


2.53 




4500 


7 


J1141-6545 


393.90 


0.20 


1.86 


0.17 


5.33 <50 


0.18 


2.30 




0.59 


8 


Unknown companion 


B1820-11 


279.83 


357.76 


200.67 


0.79 


< 10-4 


0.07 








9 


J1829+2456 


41.00 


1.17 


7.24 


0.14 


0.28 -0.02* 


0.29 


2.53 




60 
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& Weisbere (1989) 6: 


Bureav et al. (2003') 






: Wolszczan C1990'); Stairs et al. ('2002') 






7: 


Thorsett et al. ( 


19931: Arzoumanian 


(1995') 



3: Anderson et al. (1990') : Anderson (1992') 
Prince et al. (1991'): Deich fc Kulkarni (1996') 
4: Nice et al. (1996); Hobbs et al. (2003) 
5: Lvne et al. (20001 
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